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Abstract 


In this work, a fully three-dimensional mathematical model for planar porous-electrode-supported (PES) solid oxide fuel cell (SOFC) has 
been constructed to simulate the steady state electrochemical characteristics and multi-species/heat transport. The variation of chemical species 
concentrations, temperature, potential, current and current density for two types of PES-SOFC developed by central research institute of electric 
power industry (CRIEPI) of Japan are studied in the co-flow pattern. In the numerical computation, the governing equations for continuity, 
momentum, mass, energy and electrical charge conservation are solved simultaneously using the finite volume method. Activation, ohmic, and 
concentration polarizations are considered as the main sources of irreversibility. The Butler-Volmer equation, Ohm’s law, and Darcy’s gas model 
with constant porosity and permeability are used to determine the polarization over-potential, respectively. The output voltages measured in 
experiments and calculated using the above models agree well. For the cell using the same material and manufacturing process, the results show the 
type-II PES-SOFC is with better performance. However, the electrolyte of type-II PES-SOFC should be with higher maximum ionic conductivity. 
Furthermore, these results will be used to evaluate the overall performance of a PES-SOFC stack, and to significantly help optimize their design 


and operation in practical applications. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Solid oxide fuel cells (SOFC) are the basis of a promising 
technology for the high-efficiency conversion of hydrocarbon 
fuels to electric energy, and currently are considered one of 
the most promising type fuel cells used as a power source for 
stationary or distributed power plants [1,2]. SOFCs operate at 
relatively high temperature, around 700-1000 °C for maintain- 
ing high oxygen-ion-conductive of solid oxide electrolyte. Such 
high temperature accelerates electrochemical reaction; there- 
fore, they do not require precious metal catalysts to enhance the 
reaction [3,4]. Among fuel cells, SOFCs have a great advantage 
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in fuel adaptability. They are less sensitive to fuel composition 
compared to other fuel cell system. Besides the reformed gas 
from external reformer, the anode feed can employ natural gas, 
LPG, methanol, or coal gasified gas through internal reform- 
ing on the anode side. Another attractive feature is SOFC also 
can be combined with other power generation systems (e.g., 
gas turbines) to achieve high overall electrical power generation 
efficiency due to the high working temperature. 

SOFCs usually consist of an interconnect structure and a 
three-layer region composed of two porous ceramic electrodes 
(anode and cathode), and a dense, gas-tight, oxygen-ion- 
conducting electrolyte. Planar and tubular geometries are the 
two popular designs for SOFCs. The technique of designing 
and manufacturing of tubular type SOFC is relatively mature 
[1,5]. However, planar type design has received much atten- 
tion recently, because it is simpler to fabricate and easier to be 
made into various shapes than other type designs. Besides, pla- 
nar type SOFC offers higher power density relative to tubular 
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Nomenclature 

A area (m2) 

dp pore diameter in porous medium (m) 

Die overall effective diffusion coefficient of species i 
(m? s7!) 

Dij binary mass diffusion coefficient of a mixture 
species i and j (m? s7!) 

Dim mass diffusion coefficient for species i in the mix- 
ture (m2 s7!) 

Dk Knudsen diffusion (m2 s7!) 

Erpg electromotive force (V) 

F Faraday’s constant: 96,487 C mol! 

AG Gibb’s free energy change (Jkg~!) 

H total enthalpy (J kg7!) 

AH enthalpy change of the reaction (J kg7!) 

i current (A) 

J current density (A m~?) 

Jo exchange current (A m~?) 

k turbulent kinetic energy (m? s~7) 

Kn Knudesen number Kn = )¢/dp 

M molecular weight (kg kmol~!) 

Ne number of electrons 

P pressure (Pa) 

Pitpp the ratios of partial pressure over the standard 
state pressure of 1.013 x 10°Pa at the elec- 
trolyte/electrode interface ie {H2, H20, O2} 

R universal gas constant (Jkg~! K7!) 

Rohm ohm resistance (Q) 

So source or sink 

AS entropy change of the reaction (Jkg~! K7!) 

T temperature (K) 

U velocity (m s7!) 

Uf fuel utilization factor 

V potential (V) 

V°, V* electrical potential at anode and cathode, respec- 
tively (V) 

x width direction of the cell 

Xi molar fraction of species i (kmol kmol7!) 

y thickness direction of the cell 

Y mass fraction of species (kg kg~!) 

Z length direction of the cell 

Greek symbols 

a electron transfer coefficient (usually taken to be 
0.5) 

Io diffusion coefficient (kg m`! s7!) 

ô thickness (m) 

E turbulent dissipation rate (m? s73) 

EP porosity of porous medium 

no over-potential incurred by the activation polariza- 
tion (V) 

n®™  over-potential incurred by the ohmic polarization 
(V) 

Kp permeability of porous medium 

À thermal conductivity (W m`! K7!) 


Àf mean free path length (m) 

v special Fuller et al. diffusion volume 
p density (kg m~?) 

pi resistivity (Q m) 

TP tortuosity of porous medium 
Subscripts 

ca cathode 

ele electrolyte 

f fuel 

8 gas 

IC inter-connector 

Ref reference 

s porous solid 

Superscripts 

a anode 

c cathode 

e electrolyte 


type SOFC, which is ascribed to low electrical resistance due to 
shorter current paths. 

In SOFC cell, the oxygen ions formed at the cathode 
migrate through the electrolyte to the anode/electrolyte inter- 
face where they react with the hydrogen and carbon monoxide 
contained in the fuel, producing water and carbon dioxide 
while releasing electrons that flow via an external circuit to the 
cathode/electrolyte interface [6]. SOFC system involves com- 
plex multi-component transport, chemical, and electrochemical 
processes, with its operating performance strongly affected 
by corresponding transport resistances and activation barriers 
[7]. Transport resistances of gas-phase species in the porous 
electrode and oxygen ions in the electrolyte, as well as the 
activation energy barriers for electrochemical reactions can 
result in various polarizations. They are represented as con- 
centration, activation, and ohmic over-potentials. These various 
polarizations are functions of both operating conditions and 
physical properties of cell’s components. The operating con- 
ditions include temperature, pressure, and fuel and oxidizer 
concentrations. The cell properties involve materials, macro- and 
micro-structures of electrolyte and composite electrodes, which 
include the porosity, tortuosity, permeability, and thickness of 
anode and cathode, the ionic conductivity and thickness of elec- 
trolyte, and the active area and activity of electrode—electrolyte 
interface. 

A goal for SOFC is to develop economically acceptable sys- 
tems. It is expected that the capital costs can be lowered through 
advancement in cell manufacturing and cell design. An increase 
in power density is a major technical contribution to further 
cut costs. To understand these various resistances, optimize 
the cell design, and enhance the cell performance, parametric 
studies for these various polarizations in the electrolyte and com- 
posite electrodes have been widely made through experiments 
[8,9], theoretical analysis [2,5,7,10,11], and numerical simula- 
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tion [1,3,12]. Experimental approaches would be very costly and 
time consuming. In addition, it would be very difficult to measure 
directly the detailed profile of chemical species, current, elec- 
tric potential, etc. in the cell. However, the three-dimensional 
numerical simulation, based on commercial computational fluid 
dynamic (CFD) code and a user defined routine for dealing with 
electrochemical processes, can provides a better understanding 
of main phenomena that govern fuel cell performance [13], for 
example, the cell internal temperature and flow composition pro- 
files, fuel and oxidant utilization, cell power output and cell 
voltage or current output. 

In this work, based on commercial CFD software (CFX4.4) 
and a routine developed by ourselves and used to calculate the 
mass diffusion and heat transfer in the porous eletrodes for 
PES-SOFC cell, the electrochemical processes and the current 
field, a fully three-dimensional mathematical model for planar 
PES-SOFC has been constructed to simulate the flow, heat and 
mass transfer, and the electrochemical performance of two types 
of PES-SOFCs. The more detailed spatial variation of chemi- 
cal species, local over-potential, electric potential and current 
density are obtained. These results will be used to evaluate 
the overall performance of the PES-SOFC stack, and to sig- 
nificantly help optimize their design and operation in practical 
applications. 


2. Computational model 


The schematic views of two types of new planar PES-SOFC 
simulated are shown in Fig. 1. The layer of electrolyte is lami- 
nated onto the diffusible porous electrode. In a practical planar 
SOFC stack, many planar cells are mounted in a container. For 
most of the planar cells mounted in serial way, each of them 
is connected by two others. Therefore, it is possible, and very 
likely, that most of the single planar SOFCs work under the same 
operating condition. This allows us to obtain results that are very 
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Fig. 1. Schematic view of two types of planar PES-SOFC: (a) type-I and (b) 
type-II. 
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useful for evaluating the performance of a cell stack, through 
analysis of heat/mass transfer and electrochemical performance 
of one single cell and its controllable area. In SOFC cell, the 
fluid consists of multiple components, oxygen and nitrogen in 
the cathode channel and hydrogen and water vapor in the anode 
channel. The heat and mass transfer are coupled with each other. 


2.1. Flow and heat/mass transfer 


The basic set of equations for flow and heat/mass transfer 
comprises the equations for conservation of mass, momen- 
tum, energy and concentration, known as the three-dimensional 
Navier-Stokes equations. The effect of turbulence is represented 
by the k—e model. In this model, wall functions are used to 
resolve the flow near wall surfaces. The following time-averaged 
equations are used for steady state in compressible gases [14]: 


V -(pUo— IgV) = So (1) 


where ¢ represents U, k and £, H and Y. 

The model commonly used for flow in porous medium is a 
generalization of the Navier-Stokes equations and Darcy’s law 
with constant porosity and permeability. Mass diffusion coef- 
ficients are required whenever species transport equations in 
multi-component flows are solved. Diffusion in porous medium 
is usually described by molecular diffusion or Knudsen diffu- 
sion. Knudsen diffusion occurs when the diameter pores in the 
porous mediums are small compared with the mean free path 
of the gas molecules, i.e. Knudsen number Kn > 1. Molecular 
diffusion occurs when the pore diameter is large compared with 
the mean free path of the gas molecules, i.e. Kn < 1. For SOFCs, 
both Knudsen and ordinary diffusion processed have to be con- 
sidered since, in general, Kn ~ 1. The mass diffusion coefficient 
for species i in the mixture is computed as 


1— Xi 
Dim = ss 
5 (X;/Dij) 
Jj+i 
The binary diffusion coefficient Dj of a mixture species i and 
j is computed using the Fuller et al. expression [15]: 


(2) 


0.01437!-75 
1/2, 1/3 1/342 
PM,/ [v}/ +077] 


Dij = (3) 


where Mj = 2(1/M; + 1/M;)"! and v can be found in literature 
[15]. 

For straight and round pores [16], Knudsen diffusion is given 
by 


D Lt 8RT (4) 
eae p mM; 


Due to the pore’s tortuous nature and pore constrictions, the 
diffusivity is corrected by tortuosity factor, tp, and porosity 
of porous medium, ep. The tortuosity factor accounts for both 
the effects of altered diffusion path length and changing cross- 
sectional area in constrictions. Zhang and Liu [17] proposed a 
new way to predict tortuosity in catalyst pellets based on fractal 
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Table 1 
Physical properties of fuel cell components 


Type I Type II 


Geometry parameters 
Length (m) 
Width (m) 
Anode thickness (m) 
Cathode thickness (m) 
Electrolyte thickness (m) 
Inter-connector thickness (m) 


Material properties 
Thermal conductivity of anode (W m~! K~!) 
Thermal conductivity of cathode (W m~! K7!) 
Thermal conductivity of electrolyte (W m7! K7!) 
Thermal conductivity of inter-connector (W m~! K7!) 
Electric resistivity of anode (Q m) 
Electric resistivity of cathode (Q m) 
Ionic resistivity of electrolyte (Q m) 
Electric resistivity of inter-connector (Q m) 
Porosity of anode 
Porosity of cathode 
Permeability of anode 
Permeability of cathode 
Tortuosity of anode 
Tortuosity of cathode 
Pore diameter of anode 
Pore diameter of cathode 


60 x 1073 
50 x 107° 
5x 1073 
5x 1073 


60 x 107° 

10 x 1073 

10 x 1073 
0.05 x 1073 
0.5 x 1073 


6.23 
9.6 
2.7 
9.6 
2.98 x 10> exp(—1392/T) [10] 
8.11 x 1075 exp(600/T) [10] 
2.94 x 1075 exp(10,350/7) [10] 
6.41 x 10-8 [3] 
0.38 
0.5 
1.0 x 10-12 
1.0 x 10-12 
1.5 
1.5 
2.0 x 1076 
2.0 x 1076 


geometries; a hard spheres model is often used for determining 
tortuosity. The accurate determination of this parameter is out 
of the scope of this paper and the measured values in Table 1 
were used instead. The overall effective diffusion coefficient is 
given by [16]: 


ep{ 1 1 i 
Die = — +— (5) 
Tp (5 Dx 


In order to take into account the heat conduction of porous 
solid, the effective thermal conductivity in the porous medium 
is computed by, as the volume average of the fluid conductivity 
and the solid conductivity: 


À = Epdg + (1 — Ep)As (6) 


where Àg is the thermal conductivity of multi-component mix- 
ture, and its computational formula and coefficient can be found 
in literature [15]. A, is isotropic constant and its measured values 
are listed in Table 1. 

In addition, there is large resistance of flow in the porous 
medium. The isotropic linear resistance tensor is proportional to 
the inverse of the permeability «. The resistance force, i.e. body 
force, is the source of momentum equation. 


2.2. Electrochemical model 


The electrochemical model calculates the current power out- 
put and the molar compositions of cathode and anode flows for 
each control volume. For clarity and comparing with the exper- 
imental results, the calculating conditions are the same as the 
experimental conditions. The fuel consists of H? and H2O on 


the anode side and the oxidizer (air) is modeled as an O2/N2 
mixture on the cathode side. Fuel and air run in a co-current 
flow. The electrochemical reaction occurring at the interface of 
cathode/electrolyte is 


O2 + 4e7 > 2077 (7) 


Then, the oxygen ion is transported through the electrolyte 
and into the active anode. The electrochemical reaction occur- 
ring at the interface of anode/electrolyte is 


Hy + O” > H2O + 2e7 (8) 


The electromotive force yielded by the oxidation of H3 is a 
local quantity, as it depends on the gas composition and tempera- 
ture, and can be determined by the well-known Nernst equation, 
as following [1]: 


Pu, TPB Pp 
In 2 O2,TPB (9) 


-AG —-AG — RT 
= + 
PH,O,TPB 


2F 2F 2F 


E tpg = 


In the electrochemical reaction proceeding, the changes of 
chemical enthalpy at temperature T, AH, entropy, AS, and 
Gibb’s free energy, AG, have the thermodynamic relationship. 


AG = AH — TAS (10) 


Concentration over-potentials appear when the mass trans- 
port hinders the electrode reaction. Concentration polarization 
becomes an important loss at high current densities and small 
fuel concentrations (Uf > 80%). The main factors that contribute 
to concentration polarization are the diffusion of gases through 
the porous media and solution dissolution of reactants and prod- 
ucts. Along with the consumption and production of the gas 


Y. Wang et al. / Journal of Power Sources 170 (2007) 101-110 105 


species in fuel and air flows, the molar fraction of the reactants 
and products from the electrochemical reactions also vary in the 
flow stream of fuel and air. In Eq. (9), the Nernst electromotive 
force obtained includes the over-potential incurred by the com- 
ponent diffusion, due to using the species concentrations at the 
three-phase boundaries, i.e. considering the mass diffusion of 
multi-component flow in porous medium. This yields the local- 
ized electromotive force, Erpg, over the electrolyte layer, which 
yields the localized ionic transfer rate through the electrolyte 
layer. 


. Erpe nt- (VS — V8) 
l = RE 


ohm 


(1) 


2.2.1. Activation polarization 

Activation polarization is controlled by the electrode kinet- 
ics at the electrode surfaces. This polarization is directly related 
to the activation barrier that must be overcome by the reacting 
species in order for the electrochemical reaction to occur. Activa- 
tion over-potential, 7°“, incurred by the activation polarization, 
reflects the kinetics of reactions and occurs at both anode and 
cathode. At high operating temperature, the electrode reaction 
is rapid and, as a result, the activation over-potential is usually 
small. However, as the operating temperature falls, it can become 
the most significant cause of voltage drop. It is often represented 
by the non-linear Butler—Volmer equation [8], which relates the 
current density drawn to the activation over-potential and, for a 
first order charge transfer controlled electrochemical reaction, 
is given by: 


j Ne Fret Nne Fy 
> = 1.0 12 
5p (« = exp (—(.0- a)" (12) 


where jo the exchange current density, given by [8]: 


: Pu, Pio 100 x 10° 
a = 5.5 x 108 (=) ( 2 ) exp | -——— ]) 03) 
J Pal Bee RT 


Po 0.25 
j =7.0 x 108 (=) ex 
PRef 


120 x 10° 


14 
RT (14) 

The non-linear Butler—Volmer equation is solved numerically 
using the Newton—Raphson method in this work. The computa- 
tional accuracy is 1.0 x 1076. 


2.2.2. Ohmic losses 

Ohmic losses are caused by resistance to conduction of ions 
(through the electrolyte) and electrons (through the electrodes 
and inter-connector) and by contact resistance between cell com- 
ponents. This voltage drop is important in all types of cells and 
is essentially linear and proportional to current density. Because 
the ionic flow in the electrolyte and the electronic flow in the 
electrodes obey Ohm’s law, ohmic losses can be expressed by: 


n®™ = Romi (15) 


The ohm resistance is calculated according to the second 
Ohm’s law: 


ô 
Rohm = Pi A (16) 


where p; is the corresponding material resistivity, calculated 
with temperature-dependent relation (Table 1). 

In this work, the current calculation in the electrolyte is one 
dimension, along with the thickness direction, due to very small 
thickness, and it is three dimensions in the anode and cathode. 
Applying Kirchhoff’s law of current, Eq. (17), the potential and 
current in each calculating volume are obtained: 


y7=00 (17) 


where ij =(V;— Vp)/R;,p, P is the present calculating volume 
and j indicates six neighborhood volume at x-, y-, and z-axis. 
The difference between the current through every layer of the 
electrode and the given current is below 1.0 x 1074. 

Since the potential difference between two inter-connectors 
is the cell terminal voltage, the potential at the outside surface 
of inter-connector in contact to the anode is assumed to be zero, 
thus the potential at the outside surface of inter-connector in 
contact with the cathode will be the terminal voltage of the fuel 
cell. Once all the local electromotive forces are obtained from 
Eq. (9), there are only two unknown condition, the total current 
flowing out from the cell and the potential at the inter-connector 
in contact with the cathode. Therefore, either must be specified 
as the initial condition. In the present work, the total current 
is specified according to the average current density and total 
reacting area of the cell. 


3. Validity of models 


The schematic view of experiment measurement of single cell 
for type-I PES-SOFC is shown in Fig. 2. In the figure, the Ref 
indicates the measuring reference point, the potentials on the 
cathode and the anode are measured at three measuring posi- 
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Fig. 2. Schematic view of experiment measurement. 
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Fig. 3. Voltage (type-I cell): (a) voltage at open circuit and (b) output voltage 
of experiment and calculation. 


tions, and the current is collected using the platinum wire. The 
experiment measurements are completed by Dr. Ito et al. (the 
researcher of materials science research laboratory of CRIEPI) 
under 1000 °C temperature. 

The voltage at open circuit and the output voltage obtained 
from experiment and calculation at the same operating condi- 
tions and structure to the experiment equipment are given in 
Fig. 3. In Fig. 3a, it can be seen that the voltage at open circuit 
increases and trends to the constant, as the increasing of fuel flow 
rate. Then, the variation of the output voltage as a function of the 
current density under 300 N ml min™! fuel and 1000 N ml min7! 
air is measured. The output voltage decreases with the increas- 
ing of the current density. The output voltage of the calculation 
agrees well with that obtained for the experiment. This proves 
the above models used in the analysis are accurate and cred- 
itable, and is suitable to be used to analyze the performance of 
PES-SOFC cell [18]. 


4. Results and discussions 


The measured physical properties of the cell components 
are listed in Table 1. The operating conditions are listed in 
Table 2. In the calculation, the fuel utilization is 80%, and the 
oxygen utilization is 10%. In the following, the more detailed 
spatial variation of temperature, chemical species, local over- 
potential, electric potential and current density are given. Due 
to most of single planar SOFC working under the same operat- 
ing condition, the outside surfaces of two inter-connectors are 
specified to the periodic boundary. This means the heat exchange 


Table 2 
Operating conditions 
Type I Type II 

Operating pressure (Pa) 1.013 x 10° 1.013 x 10° 
Inlet temperature (°C) 1000 1000 
Inlet mass flow rate of air (N ml min!) 4775 4775 
Inlet mass flow rate of fuel (N ml min!) 130.6 108.8 
Average current density (A m~*) 4000 4000 
Electrochemical reaction area (m7) 3.6 x 1073 3.0 x 10-3 
Fuel utilization (%) 80 80 
Species molar fraction of fuel 

H2 0.9578 0.9578 

H20 0.0422 0.0422 
Species molar fraction of air 

O2 0.198 0.198 

N2 0.802 0.802 


between the cells is the same. In the calculation, there is the no- 
electrochemical-reaction extend region (4mm) at the inlet and 
outlet of cell. In the following figures, the midsection of pipe 
indicates the center-section of pipe on the porous electrodes, and 
the porous region indicates the cut-section between two pipes in 
the porous electrodes. 


4.1. Distribution of temperature 


Fig. 4 shows the temperature field for two types of PES- 
SOFC. On the whole, the fuel and air flow is progressively heated 
up by the heat generation due to the entropy change of the elec- 
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Fig. 4. Distribution of temperature: (a) type-I and (b) type-II. 
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Fig. 5. Distribution of oxygen molar fraction in the cathode: (a) type-I and (b) 
type-II. 


trochemical reaction at the interface of anode and electrolyte. 
The temperature of air flow in cathode channel is lower, espe- 
cially in the pipe of cathode, due to its larger mass flow rate. 
In the anode channel, the fuel mass flow rate is very small, so 
its velocity is very low, especially in the porous medium due 
to large resistance force. In addition, the solid thermal con- 
ductivity of porous medium relative to the fuel flow is large. 
Therefore, the temperature distribution in the porous medium of 
anode channel is mainly determined by heat conduction. Due to 
the different thermal conductivity of different material of the cell 
components, there are larger temperature drops on the interface 
between different cell components. 

From this flow arrangement, the hotspot temperature of the 
cell occurs in the output end of anode near the electrolyte. Nev- 
ertheless, the physical properties of the cell components, heat 
generation and input temperature and mass flow rate of air and 
fuel will collectively affect the temperature field in the fuel cell. 
Thus the hotspot position and temperature will depend on the 
operating condition of the fuel cell. 


4.2. Distribution of species 


Fig. 5 shows the distribution of oxygen molar fraction in 
the cathode channel. It is observed that the oxygen molar frac- 
tion decreases along the air flow direction due to the oxygen 
reduction reaction. The contour shape of oxygen also indicates 
a relatively large difference of oxygen molar fraction between 
the bulk flow and the cathode/electrolyte interface. Comparing 
Fig. 5a with b, the local molar fraction of oxygen in the cath- 
ode channel of type-II PES-SOFC is higher, due to the smaller 
distance between the bulk flow and the cathode/electrolyte inter- 
face. Fortunately, the mass flow rate of air is enough in order to 
control the cell temperature and prevent the cell broken. Thus, 
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Fig. 6. Distribution of hydrogen molar fraction in the anode: (a) type-I and (b) 
type-II. 


the mass transport resistance in the cathode channel might have 
an unimportant effect on the cell performance. 

The distributions of hydrogen and water vapor molar fraction 
in the anode channel are shown in Figs. 6 and 7, respectively. 
Since the hydrogen consumption occurs at the interface of 
anode/electrolyte, the hydrogen molar fraction decreases along 
the fuel stream. Corresponding to this situation for hydrogen, the 
production due to the electrochemical reaction makes the water 
vapor increases gradually along the fuel stream. The shape of the 
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Fig. 7. Distribution of water vapor molar fraction in the anode: (a) type-I and 
(b) type-II. 
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contour lines of hydrogen and water vapor in the fuel channel 
is firstly very sharp and then relatively flat from the interface 
of anode/electrolyte to the interface of anode/inter-connector, 
and this also implies that species diffusion over-potential at the 
anode is larger than it at the cathode. Comparing Fig. 6a with b 
and Fig. 7a with b, the shape of the contour lines of hydrogen and 
water vapor in the fuel channel of type-II PES-SOFC is flatter 
relative to type-I PES-SOFC. This is owed to the smaller distance 
between the interface of anode/electrolyte and the pipe of porous 
medium. In a word, the mass transport resistance on the fuel 
channel might be dominant in lowering the cell performance, 
especially for high fuel utilization of the cell. 


4.3. Distribution of current field 


Fig. 8a and b shows the electric potential distributions on 
the mid-planes of the type-I and type-II PES-SOFC along the 
flow direction, respectively. In the porous medium of anode, 
due to the outside surface of inter-connector in contact to the 
anode assumed to be zero, the electric potential is the negative 
value, and it reduces from the outside surface of inter-connector 
to the interface of anode/electrolyte. Furthermore, due to the 
high electric conductivity of the anode material, its absolute 
values is small, the absolute maximum at the interface of 
anode/electrolyte, i.e. the ohmic over-potential of anode and 
inter-connector in contact to the anode is small. In the porous 
medium of cathode, the electric potential is the positive value, 
and it reduces from the interface of cathode/electrolyte to the 
outside surface of inter-connector in contact to the cathode. The 
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Fig. 8. Distribution of electric potential: (a) type-I and (b) type-II. 


electric potential at the outside surface of inter-connector in con- 
tact to the cathode is the terminal voltage. At the region between 
two pipes of porous electrodes, the current through that area 
is relatively large which leads to more dense contour lines, i.e. 
there is a large drop of electric potential. It is also seen that the 
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Fig. 9. Distribution of Nernst potential: (a) type-I and (b) type-II. 
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Fig. 10. Distribution of current density at the electrolyte: (a) type-I and (b) type-II. 


cathode ohmic over-potential, the difference of electric potential 
between the interface of cathode/electrolyte and the outside sur- 
face of inter-connector in contact to the cathode, is larger than 
the anode ohmic over-potential. Comparing Fig. 8a with b, the 
ohmic over-potential of cathode and anode of type-II PES-SOFC 
is lesser than that of type-I PES-SOFC, due to the smaller thick- 
ness, and the terminal voltage of type-II PES-SOFC is larger 
than that of type-I PES-SOFC. 

Fig. 9 shows the distribution of electromotive force, i.e. 
Nernst potential, taking into account the diffusion polarization. 
Due to the progressive decrease of hydrogen and oxygen concen- 
tration, the Nernst potential decreases along the flow direction 
(z-axis). The different diffusion of species in the pipe and porous 
region of the anode yields the Nernst potential presents the wave 
variation along x-axis. On the whole, the Nernst potential of 
type-II PES-SOFC is larger than that of type-I PES-SOFC, due 
to its smaller thickness of anode. 

Fig. 10 shows the distribution of current (ionic transfer rate) 
density through the electrolyte. According to Eq. (11), the cur- 
rent density distribution is determined by many parameters, for 
example, temperature, species concentrations, potential of anode 
and cathode, activation over-potential, etc. As the temperature 
increases from inlet to outlet of the cell due to the electrochem- 
ical reaction, the ionic resistance of the electrolyte decreases. 
In addition, due to the larger heat exchange in the pipes of 
porous cathode, there is lower temperature at the regions of 
the electrolyte corresponding to these pipes, i.e. there is higher 
ionic resistance at these areas. Therefore, the ionic resistance 


of the electrolyte also presents the wave variation along x-axis, 
especially, for the type-I PES-SOFC. Because there is no the 
electronic conductor in the hollow pipes of porous electrodes, 
the potential and current in these pipes region is zero. The cur- 
rent accumulated at the interface of cathode/electrolyte moves 
around the hollow pipes region to the inter-connector. Con- 
versely, the current in the anode moves form the inter-connector 
in contact to the anode to the interface of anode/electrolyte. The 
potential V° of cathode near the electrolyte between these pipes 
and the electrolyte is higher in the cathode, whereas V° is lower 
in the anode. Therefore, the current density also waves along 
the x-axis. There is lower value at the region corresponding to 
these pipes. Comparing Fig. 10a with b, the current density at 
the electrolyte of type-I PES-SOFC is evener than that of type-II 
PES-SOFC. 


5. Conclusions 


In the work, a fully three-dimensional numerical analysis of 
the fluid dynamics and electrochemical kinetics is conducted 
for two types of planar PES-SOFCs. Significant results about 
the spatial variation of temperature, species concentration, elec- 
tric potential, etc. are presented. It is found that the fuel and air 
flow is progressively heated up along with the flow direction, 
and the hotspot locates at the end of the anode near the elec- 
trolyte, when the air and fuel stream is the co-flow. The mass 
transport resistance of species in the porous anode of type-I PES- 
SOFC under higher fuel utilization is dominant in lowering the 
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cell performance due to the larger thickness of anode. Under 
the same fuel utilization and average current density, the out- 
put electric potential of type-II PES-SOFC is higher than that 
of type-I PES-SOFC. This also means the cell performance of 
type-II PES-SOFC is better than that of type-I PES-SOFC. How- 
ever, the electrolyte of type-I] PES-SOFC should be with higher 
maximum ionic conductivity. 
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